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1.  Introduction 

A  non-perturbative  general  relativistic  approach  to  global  astrometry  was  de¬ 
veloped  by  de  Felice  et  al.  (1998)  to  handle  satellite  astrometry  data  in  a  gen¬ 
uine  relativistic  framework.  In  this  contribution,  the  framework  above  has  been 
further  exploited  to  account  for  stellar  motions  and  parallax.  Because  of  the 
relevance  that  accurate  knowledge  (to  10-5  or  better)  of  the  relativistic  parame¬ 
ter  7  has  to  fundamental  physics,  a  Parametrized  Post-Newtonian  (PPN)  model 
has  also  been  implemented,  which  allows  the  direct  estimation  of  7  along  with 
the  astrometric  parameters.  These  models  have  been  tested  on  end-to-end  sim¬ 
ulations  of  the  mission  GAIA.  The  results  show  that,  within  the  limitation  of 
the  simulation  and  the  assumptions  of  the  adopted  model,  measurements  accu¬ 
rate  to  100  juarcsec  of  large  arcs  among  stars  repeated  over  a  few  years  can  be 
modelled  to  establish  a  dense  reference  frame  with  a  precision  of  a  few  tens  of 
juarcseconds.  Moreover,  our  experiments  indicate  that  7  can  be  estimated  to 
better  than  10-6. 


2.  Relativistic  effects  in  GAIA 

In  the  data  reduction  of  the  Hipparcos  astrometric  mission,  the  observations 
were  pre-corrected  for  relativistic  effects  to  (w/c)2  —  GMq/c2Rq  ~  2  mas.  The 
observation  equations  were  then  formulated  in  a  three-dimensional  Euclidean 
space.  The  expected  astrometric  precision  for  the  GAIA  mission  is  of  the  order 
of  100  juarcsec  for  stars  of  magnitude  V  =  17,  and  ~  10  ^arcsec  for  V  =  12 
(Lindegren  and  Perryman,  1996).  To  meet  such  goal,  it  is  mandatory  that  the 
observations  be  modelled  at  the  level  of  ~  1  ^arcsec,  in  a  framework  which 
utilizes  a  relativistic  theory  of  gravitation. 

Especially  significant  for  the  GAIA  mission  are  metric  perturbations  due  to 
the  solar  system  planets  which  will  affect  the  path  of  the  photons  reaching  the 
satellite.  The  following  table  lists  the  light  deflection  angles  due  to  the  major 
planets  at  their  closest  approach  to  Earth,  and  for  different  elongations  from  the 
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light  source.  We  have  not  considered  the  internal  planets  Venus  and  Mercury 
because  they  would  be  too  close  to  the  Sun  to  be  observed  by  GAIA.  The  mass 
of  the  planet  is  denoted  by  iff,  r®  is  the  minimum  distance  planet-Earth,  D  the 
corresponding  angular  diameter  of  the  planet  and  Sa  the  light  deflection  angle. 
It  can  be  seen  that  the  mass  of  Jupiter  produces  significant  effects  up  to  angular 
distances  of  ~  15°,  while  Saturn  needs  to  be  considered  up  to  ~  2°.  The  mass 
of  Uranus  and  Neptune  might  be  important  only  for  the  brightest  stars. 


Table  1.  Deflections  of  a  light  ray  due  to  solar  system  planets  as  seen 
from  the  Earth  for  a  limb-grazing  ray  ( Saig ),  and  for  2  and  15  degrees 
of  elongation. 


Planet 

M  •  10_M 
kg 

r®  •  10“b 
km 

D 

(") 

$  C^lg 

mas 

Sa( 2°)  |  Sa(  15°) 

fi&Tcsec 

Mars 

6.4191 

78.3 

8.94 

0.116 

— 

— 

Jupiter 

18992 

628.7 

23.43 

16.289 

53.0 

7.0 

Saturn 

5686.5 

1277.4 

9.76 

5.765 

7.8 

— 

Uranus 

868.49 

2720.0 

1.80 

2.235 

— 

— 

Neptune 

102.35 

4347.0 

1.06 

2.800 

- 

- 

3.  Modelling  the  observations 

The  adopted  mathematical  model  consists  of  a  satellite  (the  observer  orbiting  on 
a  spatially  circular  geodesic  around  the  Sun.  The  latter  is  assumed  non  rotating, 
spherically  symmetric  and  therefore  generating  a  Schwarzschild  space-time  met¬ 
ric,  the  Sun  conciding  with  the  barycenter  of  the  satellite-Sun  system.  In  this 
context  it  is  convenient  to  adopt  polar  coordinates:  the  colatitude  0  £  [0,7r],  the 
azimuth  0  £  [0,27t]  and  the  polar  distance  r  £]0,  oo[.  When  the  Schwarzschild 
metric  is  expressed  in  the  Parametrized  Post-Newtonian  (PPN)  formalism  and 
isotropic  coordinates,  the  two  new  parameters  7  and  (3  are  explicitely  intro¬ 
duced  (Misner  et  al. ,  1973).  These  parameters  are  equal  to  unity  in  the  theory 
of  General  Relativity. 

The  observable  quantity  for  a  GAIA-like  mission  is  the  light  coming  from 
two  objects  in  different  directions  across  the  sky  reaching  the  satellite  at  the 
same  proper  time  r.  In  particular,  the  cosine  0i2  of  the  angle  between  a  star 
pair  is  given  by  (Brumberg,  1991): 


COS  -012  = 


hapk^k^ 

y/hl7TkL1k^Jhpak%k?, 


(1) 


where  k\  and  k2  are  the  null-geodesics  from  the  two  stars,  hap  =  gap  +  uaUf3 ,  is 
a  tensor  projecting  in  the  rest  frame  of  the  observer,  and  ua  is  the  four- velocity 
vector  representing  the  observer’s  trajectory.  Our  aim  is  to  express  the  observ¬ 
able  as  function  of  same  space  coordinates,  as  well  as  of  the  PPN  parameters 
in  the  post-Newtonian  approximation.  Equation  1  does  not  explicitely  contain 
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such  parameters,  and  they  are  introduced  by  integrating  the  path  of  the  photons 
from  the  light  source  to  the  observer  (de  Felice  et  al. ,  2000).  Finally,  we  obtain 
an  equation  of  the  form: 

cos  ^12 (r)  =  /  (ri (/),  01  (/),  0r(/),  r2(Z),  02(Z),  02(Z);  7,  /?)  •  (2) 

We  note  that  the  time  at  which  the  angular  separation  of  the  star  pair  is 
observed  is  the  proper  time  r  of  the  observer.  Since  the  catalog  values  of  the 
star’s  parameters  are  referred  to  coordinate  time,  we  have  to  convert  r  into  t  in 
order  to  compute  the  correct  stellar  position. 

The  polar  coordinates  (r,  0 ,  <jr)  (the  only  ones  appearing  as  arguments  of 
the  function  /  in  a  non-perturbative  approach)  are  functionally  related  to  the 
classical  quantities  parallax  and  proper  motions,  as  described  in  the  following 
section. 


4.  The  condition  equation 

The  motion  of  a  distant  star  at  coordinate  time  t  can  be  formally  described  by 
using  Taylor  expansions  with  respect  to  an  intial  time  to  as: 

r(t)  =  r(t0)  +  r(t  -  t0)  +  0.5  r(t  -  t0 )2  +  ... 

<  9{t)  =  9  (to)  +  9(t  —  to)  +  0.5  9(t  —  to)2  +  ...  (3) 

4>(t)  =  <f>(t0)  +  4>(t  —  to)  +  0.5  <j)(t  —  to)2  +  ... 

In  principle,  r,  9 ,  and  cj)  are  not  directly  observable;  however,  since  r  >>  r®, 
the  above  coordinates  can  be  used  to  define  the  observable  stellar  parallax  and 
proper  motion ,  i.e.,  p  =  r®/r,  r®  =  1 AU ,  pe  =  9,  =  <f>.  To  conveniently 

truncate  the  Taylor  expansions,  let  us  evaluate  the  rate  of  change  of  p,  p ^  and 
pg,  i.e.,  p  =  —  pg  =  9 ,  and  p ^  =  <f>.  It  can  be  seen  that  (Green,  1985)  p  = 
—  0.2  Vrp2  sinl",  and  p  =  —0.1  VrVtp 2  sin  1",  where  p  =  yjpg  +  p ^  is  expressed 
in  arcseconds/yr,  Vr  and  Vj  are  the  star’s  bary centric  radial  and  tangential 
velocities  in  km/sec,  and  p  is  in  arcseconds.  Straightforward  calculations  show 
that,  while  p  is  negligible  in  virtually  all  cases,  the  effect  of  p  will  be  detected  by 
GAIA  for  high-velocity  stars  (~  100  Km/sec)  in  the  vicinity  of  the  Sun  (~  100 
pc),  which  represent  a  very  low  fraction  of  the  total  number  of  potential  targets. 
Therefore,  for  the  purpose  of  this  experiment,  the  final  model  included  the  zero- 
th-order  term  in  r  and  up  to  the  first-order  terms  in  9  and  cj)  of  Equation  3.  To 
estimate  the  astrometric  parameters,  along  with  the  PPN  parameters,  from  the 
measured  arcs,  we  substitute  Equation  3  in  2,  linearize  Equation  2  with  respect 
to  some  a  priori  known  values,  then  solve  the  linear  system  using  a  least-squares 
method.  The  resulting  condition  equation  is  the  following: 

2 

-  sin  ^12^12  =  +  Ci6(f>i  +  Dl8pgl  +  E^p^)  +  FSj  +  G8/3 , 

8  =  1 

where  the  coefficients  A®  C),  Tk,  A®  F  and  G  are  partial  derivatives  of 
the  function  /  calculated  at  catalog  values.  The  quantity  <5^12  is  the  difference 
between  the  measured  and  the  catalog  value  of  the  angular  separation  between 
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the  star  pair.  The  small  differences  in  the  right-hand  side  of  the  condition 
equation  represent  corrections  to  the  sought  for  parameters. 

At  present,  we  have  only  developed  the  calculations  to  estimate  8j,  though 
the  inclusion  of  80  does  not  represent  a  major  difficulty  in  principle. 


5.  Simulation  of  the  observations 

The  end-to-end  simulation  of  the  GAIA  observations  follows  the  same  three- 
step  procedure  used  for  the  static  case  (de  Felice  et  al. ,  1998),  with  the  obvious 
complications  due  to  the  inclusion  of  stellar  proper  motions  and  parallaxes,  and 
the  introduction  of  the  relativistic  parameter  7.  Measurements  (arcs  between 
stars)  are  generated  by  a  satellite  that  sweeps  the  sky  according  to  a  Hipparcos- 
like  scanning  law,  and  the  option  of  three  viewing  directions  (FOVs),  studied 
for  GAIA,  is  implemented. 

Stars  are  simulated  within  a  uniform  density  sphere  of  radius  500  pc  (i.e. 
p  >  2  mas)  and  proper  motions  are  computed  for  each  star  on  the  basis  of  its 
radial  velocities,  randomly  generated  with  (ty)  =  0  km/s  and  oVr  =  15  km/s. 
Table  2  lists  the  numerical  values  adopted  for  the  main  simulation  parameters. 


Table  2.  Relevant  parameters  of  the  GAIA  dynamic  simulation. 


parameter 

numerical  value 

satellite  orbital  radius 

1.496-  1011  m 

satellite  precession  angle 

O 

CO 

satellite  spin  period 

128  min 

angles  between  FOVs 

54°,  78°5 

amplitude  of  each  FOV 

1°6 

radius  of  the  simulated  sphere 

2  mas 

mean  stars  transverse  velocity 

0  km/s 

0  of  the  transverse  velocity  distribution 

15  km/s 

catalog  error  for  parallaxes 

2  mas 

catalog  error  for  angular  coordinates 

2  mas 

catalog  error  for  proper  motions 

2  mas 

6.  Results  and  discussion 

Table  3  shows  the  final  astrometric  errors  for  seven  different  simulations  of  the 
non-perturbative  dynamical  model.  The  mission  length  varies  from  one  to  five 
years.  The  catalog  reference  (mean)  time  is  always  at  half  the  mission  lifetime 
(to  =  AT/2),  which  provides  the  most  accurate  positions,  while  the  correlation 
between  position  and  proper  motion  is  minimized  (Eichhorn  and  Googe,  1969). 
We  note  that  the  true  errors  of  the  five  stellar  parameters  scale  well  with  the 
square  root  of  the  number  of  observations  (n0&s).  It  is  also  evident  that  the  model 
becomes  inadequate  below  AT  =  2  yr,  as  the  relatively  short  time  baseline  makes 
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it  difficult  to  disentangle  the  proper  motion  components  of  stellar  motion  from 
the  parallax. 

Table  4  reports  the  astrometric  results  along  with  the  estimation  of  7  in  the 
perturbative  (PPN)  case  for  a  static  sphere  (only  star  positions  are  considered), 
and  a  mission  duration  of  1  year.  With  these  runs  we  have  tested  the  influence 
of  the  number  of  stars  (JV*)  on  the  estimation  of  7.  The  results  show  that  with 
a  measurement  error  of  10  ^arcsec,  an  increase  of  IV*  from  2000  to  5000  gives 
6p  in  the  range  5  •  10-5  —  6  •  10-6.  These  results  made  us  confident  that  even 
better  results  will  be  achieved  by  GAIA;  the  inclusion  of  the  PPN  parameter  (3 
can  be  done  under  the  same  framework  with  a  little  increase  of  mathematical 
complexity.  Further  experiments  in  this  direction  are  underway. 


Table  3.  Astrometric  errors  for  the  non-perturbative  dynamic  case. 
Additional  input  parameters  are:  aoss  (single-measurement  error)  = 
100  n arcsec  and  IV*  (number  of  simulated  stars)  =  2000. 


AT 

(yr) 

H0bs 

® Sp 

(pi  as) 

&S0 

(pia,s) 

^"sin  6  6(f) 
(piSuS) 

a Sue 

(/Wyr) 

^"sin  6  6/t^ 

(/Wyr) 

5 

286431 

15.79 

9.49 

13.54 

6.58 

7.84 

4 

227869 

17.88 

10.67 

15.28 

9.60 

11.31 

3 

172150 

20.73 

12.97 

17.88 

15.65 

19.34 

2.5 

144063 

22.98 

15.56 

19.91 

21.52 

25.32 

2 

115218 

29.75 

21.06 

26.64 

67.32 

52.69 

1.5 

86028 

262.92 

200.75 

136.83 

679.07 

751.66 

1 

57892 

3691.33 

1154.09 

1927.44 

5930.67 

6827.99 

Table  4.  Final  errors  of  astrometric  parameters  and  the  PPN-7  pa¬ 
rameter  for  the  perturbative  (PPN)  static  case.  Additional  input  pa¬ 
rameters  are:  aoss  =  10  ^arcsec  and  AT  =  1  yr. 


A* 

H0bs 

(60) 

(/Wyr) 

&80 

(/Wyr) 

(sin  9  6(f)) 

(/Wyr) 

^"sin  6  6(f) 

(/Wyr) 

6p  •  10a 
(it  ~  7c) 

2000 

58987 

-0.08 

4.69 

-1.52 

10.95 

0.0532 

2500 

92061 

-0.01 

2.09 

-1.87 

2.72 

0.0090 

3000 

132402 

-0.01 

2.37 

-2.08 

2.51 

0.0364 

3500 

177642 

0.06 

2.40 

2.26 

2.34 

-0.0269 

4000 

230039 

0.25 

1.59 

0.59 

1.90 

-0.0071 

4500 

289399 

0.03 

1.63 

-2.03 

1.78 

-0.0097 

5000 

358002 

-0.07 

2.95 

3.90 

2.46 

-0.0066 
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